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Abstract 

While individual water molecules adsorb strongly on a talc surface (hydrophilic behavior), a 
droplet of water beads up on the same surface (hydrophobic behavior). To rationalize this di- 
chotomy, we investigate the influence of the microscopic structure of the surface and the strength 
of adhesive (surface-water) interactions on surface hydrophobicity. We show that at low relative 
humidity, the competition between adhesion and the favorable entropy of being in the vapor phase 
determines the surface coverage. However, at saturation, it is the competition between adhesion 
and cohesion (water-water interactions) that determines surface hydrophobicity. The adhesive 
interactions in talc are strong enough to overcome the unfavorable entropy, and water adsorbs 
strongly on talc surfaces. However, they are too weak to overcome the cohesive interactions, and 
water thus beads up on talc surfaces. Surprisingly, even (talc-like) surfaces that are highly adhe- 
sive, do not fully wet at saturation. Instead, a water droplet forms on top of a strongly adsorbed 
monolayer of water. Our results imply that the interior of hydrophobic zeolites suspended in water 
may contain adsorbed water molecules at pressures much smaller than the intrusion pressure. 
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A. Introduction 

Wetting properties of minerals in soils and rocks play a crucial role in the transport, and 
thus availability, of water and oil. Clay minerals are particularly interesting, not only due 
to their abundance in nature and in synthetic materials, but also because the existence of 
clays with different structures allows us to investigate the effect of surface microstructure 
on macroscopic properties. Clay surfaces can be either charge-neutral or have a net charge, 
which is balanced by counter-ions in solution. Molecular simulation has furthered our un- 
derstanding of both these types of clays: uncharged clays have been studied using both 
ab-initio^ ^ and classical simulations^, whereas simulations of charged clays have provided 
insights into interlayer propertieS^, swellin^SHEl^ g^j^d cation exchang#^^. These studies 
have shown that the surface microstructure is expected to be more important in determining 
surface-water interactions in uncharged clay d^^ ' ^^ ', and it is these surfaces that are the focus 
of the current work. Among uncharged clays, talc surfaces have attracted a lot of atten- 
tioii^^^, because of their peculiar behavior with respect to water. Water adsorption at low 
relative humidity (RH) reveals the presence of strong binding sites on talc^^. Such strong 
binding sites are absent in other uncharged clays such as pyrophyllite and fluorotalc. Yet, 
experimental contact angles indicate that the surface of talc monocrystals is hydrophobic, 
similar to that of pyrophyllitS^IESI 

To investigate this dichotomy, here we employ molecular dynamics simulations combined 
with recently developed algorithmP^^^. In agreement with experiments, we find that at low 
RH, talc surfaces display hydrophilic behavior as water adsorbs strongly to the binding sites 
on the surface. However, at saturation, cohesive interactions dominate and the interaction 
between the surface binding sites and water is minimal, resulting in a hydrophobic surface. 

To further explore the role of surface microstructure and the strength of the adhesive 
interactions on surface hydrophobicity, we also study similar clay minerals, pyrophyllite and 
fluorotalc, as well as modified talc surfaces with a range of binding site polarities, both at 
low relative humidity and at saturation. We find that the dual hydrophilic-hydrophobic 
behavior observed in talc, is generically expected to manifest for surfaces whose adhesive 
interaction energy lies in a special range. If the adhesion to water is strong enough to 
overcome the entropy of being in the vapor phase at low RH, water adsorbs strongly to 
the surface (hydrophilic behavior). At the same time, if adhesion is too weak to overcome 
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the cohesive interactions in water, the surface is hydrophobic at saturation. For modified 
talc surfaces with strong enough adhesion to overcome the cohesive interactions, aU surface 
binding sites are occupied by water molecules at saturation, as expected. Surprisingly, 
instead of observing complete wetting, we find that a water droplet sits atop the adsorbed 
water monolayer. 

B. Microscopic Models 

Talc, fiuorotalc and pyrophyllite are uncharged clay minerals, i.e., layered silicates of 
magnesium (Mg) or aluminum (Al). They belong to the family of TOT clays: each clay 
sheet consists of a layer of octahedrally coordinated Mg or Al oxide between two layers 
of tetrahedral silicon oxide (see [l](a) - side view). The surface of these sheets displays 
hexagonal rings of Si02 tetrahedra. In talc and fiuorotalc, all octahedral sites are occupied 
by Mg atoms, while in pyrophyllite two third of these sites are occupied by Al atoms (see[l](a) 
- top view). The charge on Mg and Al is balanced by hydroxyl groups in the center of the 
hexagonal cavities. In talc, these hydroxyl groups are oriented perpendicular to the surface, 
and can participate in hydrogen bonds with water. In pyrophyllite, the hydroxyl groups are 
oriented parallel to the surface, and in fiuorotalc, they are replaced by fiuorine atoms. The 
atomic coordinates for the unit cells of these clays have been included as Supplementary 
Information. 

We use the CLAYFF force fielcP to model the interactions of the clay atoms and the 
SPC/E model to describe wateil^. Lorentz-Berthelot combination rules are used to deter- 
mine the pair Lennard- Jones parameters and a rigid clay structure is assumed. As there are 
no parameters for fiuorine in CLAYFF, we assigned it a charge equal to that of the -OH 
group in talc (-0.525) and Lennard- Jones parameters of the fiuoride ion reported in Ref.'^H. 
All simulations were performed in the NVT ensemble using the LAMMPS simulation pack- 
ag^^ at a temperature, T = 300 K, maintained using a Nose-Hoover thermostat!^. SHAKE 
was used to integrate the motion of the rigid water molecules^^ and long-range electrostatic 
interactions were computed using Ewald summation. 
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FIG. 1: (a) Microscopic clay structure (Red: O, White: H, Yellow: Si, Green: Al, Cyan: 
Mg atoms). The side and top views of the pyrophyllite clay sheet show the hydroxyl (-0H) 

groups that are parallel to the sheet. In talc (top view shown), the -OH groups are 
perpendicular to the sheet and can participate in hydrogen bonds with water. In fluorotalc 
(not shown), the talc -OH groups are replaced by F atoms, (b) Part of the simulation 
setup for studying the clay - water interface. The blue box is the observation volume, v, 
used to probe density fluctuations, (c) Simulation setup for determining contact angles. 

C. Methods 

1. Clay - water interface 

A clay-water interface is representative of the situation at saturation. The setup shown 
in [l](b) is used to calculate the local water density, p(^), as well as the water density 
fluctuations near the interface. The potential of mean force, J-", for bringing a water 
molecule from bulk to a distance z from the plane of the Mg atoms for talc and fluorotalc 
(and Al for pyrophyllite) is related to p{z) by J^{z) = — /cBTln[p(z)/pb], where /cb is the 
Boltzmann constant and pb is the bulk water density. To quantify density fluctuations, we 
measure the probability distribution, Py{N)^ of flnding water molecules in an observation 
volume V, adjacent to the clay surface, using the indirect umbrella sampling (INDUS) 
methocP^'^H w/q chose a rectangular parallelopiped of dimensions 15xl5x3A^ placed 
near the surface [seeflFb)], as the observation volume. The exact z-position of v was chosen 
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so that the mean water density in v is equal to pb- The simulation box also contained a 
fixed wall of repulsive WCA particles (not shown), placed at the top of the box (far from 
v) to nucleate a vapor-liquid buffering interface. 

2. Contact angle 

The simulation setup for contact angle measurements is shown in [l](c). The contact 
angle is determined by computing water density maps in the plane of the center-of-mass of 
the drop. The curve with density equal to half of the bulk density is then fit to a circle and 
the angle between the tangent to this circle at = 7 A and the horizontal axis is taken to 
be the contact angle. While the exact quantitative value of the contact angle depends on 
the choice of Zs^ our qualitative findings do not. 

5. Water vapor adsorption 

The adsorption of water vapor at low RH corresponds the interaction of an isolated 
water molecule with the surface. To determine the corresponding adsorption free energy, 
^Mads5 we compute J-'{z) using umbrella sampling, with the weighted histogram analysis 
method (WHAM) being used to reconstruct J-'{z) from the biased trajectories. 

D. Hydrophobicity at low and high RH 

Using the various molecular measures of hydrophobicity described above, we study talc, 
as well as fiuorotalc and pyrophyllite surfaces, both at saturation and at low RH. 

1. High RH 

Theor jPStM! g^j^^j simulation d^^ * ^^^ have shown that the mean water density near a surface 
is not a good measure of its hydrophobicity. Instead, fiuctuations away from the mean, and 
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in particular, the rare fluctuations^^ indicating the cost of creating a cavity at the interface 
correlate quantitatively with the contact angld^. Patel et al. have shown that hydrophobic 
surfaces display an enhanced probability of density depletion or a low fat tail in the Pv{N) 
distribution, while Pv{N) near hydrophilic surfaces is similar to that in bulk wateil^S. As 
shown in|2|(a), Py{N) near all three clay surfaces displays a low fat tail, indicating that 
these surfaces are hydrophobic. A slight lifting of the fat tail from talc to fluorotalc and 
pyrophyllite suggests a corresponding marginal increase in hydrophobicity. 
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FIG. 2: (a) The probability, Py{N)^ of observing A^ water molecules in a probe volume 
(v = 15 X 15 X 3 A^) displays a low A^ fat tail when v is near the surface of talc (black), 

fluorotalc (red), and pyrophyllite (blue), as compared to that when v is in bulk water 
(green), (b) Water droplet proflles corresponding to p(r, z) = 0.5pb are shown for the clay 
surfaces. The contact angles for the surfaces are similar: 96° for talc, 103° for fluorotalc, 
and 105° for pyrophyllite (based on tangents drawn at zs = 7A). (c) Potential of mean 
force, J-'(^), for the adsorption of an isolated water molecule (low RH) to the clay surfaces. 
The hydrogen atoms of the talc -OH groups are located at z = 2 A and can participate in 
hydrogen bonds with water molecules, (d) J-'{z) at the clay - liquid water interface 
(saturation). To maximize H-bonding with other waters, the binding site is no longer 

occupied. 

Another way to probe surface hydrophobicity is by simulating a sufliciently large water 
droplet on the surface and estimating the corresponding contact angle. [2](b) shows the 
average shape of droplets on the clay surfaces. The curve corresponding to p(r, z) = 0.5pb 
is a circle in the (r, z) plane, where r is the distance from the axis that passes through the 
center of mass of the droplet. The contact angles obtained by tangents drawn at zs = 7 A on 
the three surfaces are similar (talc: 96°, fluorotalc: 103° pyrophyllite: 105°), and clearly 
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indicate hydrophobic behavior. 

Rehable experimental estimates of the contact angle of water droplets on both talc 
and pyrophyllite monocrystals are between 80° and 85°* ^^ ' ^^ '. The reported values for 
measurements on powders are usually smaller due to the presence of hydrophilic sites on 
the edges of finite clay particleJ^. To the best of our knowledge, no experimental contact 
angles have been reported for fiuorotalc. For both talc and pyrophyllite, the contact angles 
obtained from our simulations (96° and 105° respectively) are somewhat larger than the 
experimental estimates, suggesting that surfaces modeled with the CLAYFF model are too 
hydrophobic. Nevertheless, amongst various commonly used clay force field^^^^, we find 
that the correspondence with experiments is closest for CLAYFF. A comparison of these 
force fields with experiments is provided in the Supplementary Information. 



2. Low RE 

To investigate the wetting behavior of clay surfaces at low RH, we calculate the potential 
of mean force, T[z\ for the adsorption of an isolated water molecule. T[z) displays a 
minimum near all the clay surfaces [see[2](c)], corresponding to an adsorption (or binding) free 
energy, A/Xads- For talc, A/Xads ^ —5.9 kcal/mol, or 10 /cbT, consistent with the formation of 
a hydrogen bond between the water molecule and the hydroxyl group in talc. In fiuorotalc, 
the hydroxyl group is replaced by fiuorine, resulting in a reduction in A/Xads to -3.5 kcal/mol. 
It also shifts the location of the minimum out by 1 A as the water is no longer strongly 
bound to the surface. Pyrophyllite, with the hydroxyl group parallel to the surface has an 
even smaller A/Xads ^ —2.8 kcal/mol, and the minimum is shifted out even more. 

To compare our estimate of A/Xads from simulations to experimental data, we ana- 
lyzed the data of Michot et al?^ using a Langmuir model. This model assumes that 
there are no interactions between the adsorbed molecules and predicts a surface coverage, 
B = (P/P*)/(l + P/P*). P* is the pressure at which half of the surface sites are occupied 
and is related to A/iads through 

where amax ^4.2 nm~^ is the surface site density, 5 ?^ 1 — 2 A is the width of the surface 
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layer, i.e. the width of the PMF weU in[2](c), and 1/(3 = k^T is the thermal energy. 

In the very low RH limit, corresponding to single water adsorption, we can safely assume 
that water molecules do not interact with each other. In this regime, B P/P* and the data 
in Figure 11 of Ref.^^, allow us to obtain an experimental estimate of P* 0.056Psat for the 
talc surface. Here, Pgat = 30 mbar is the saturation pressure of water. Using this value of P* 
in equation [l| we get an experimental estimate of A/Xads ^ ~8 kcal/mol^^. This somewhat 
stronger adsorption than that predicted from simulations using CLAYFF (-5.9 kcal/mol), is 
consistent with the overestimate of the CLAYFF talc contact angle. 

If we further assume that the adsorbed water molecules do not interact with each other 
even at higher RH, the Langmuir model (with P* = 0.056Psat) predicts that B 0.9 at 50% 
RH. As water coverage on the talc surface can be large even at moderate RH, interactions 
between water molecules may be important, consistent with suggestions that clustering needs 
to be considerecF^^. In contrast, for fluorotalc B at saturation estimated from A/Xads is very 
small (?^ 1.5%), in agreement with the hydrophobic adsorption behavior reported in Figure 
10 of Ref.lSl. 

Thus, the clay surfaces simulated using the CLAYFF force field are more hydrophobic 
than the real clay surfaces used in experiments. However, the interesting dichotomy of talc 
surfaces is also observed in the simulations and our findings are qualitatively consistent with 
the experiments, both at low RH (strong adsorption for talc and not the other clays) and 
at high RH (large contact angles for all clays). 

E. Cohesion vs Adhesion 

To investigate the disparate behavior of talc surfaces at low and high RH, we compare 
T{z) for moving a water molecule away from the surface under both conditions. At satu- 
ration, J-{z) for the clay surfaces are similar (2](c)], consistent with similar droplet contact 
angle on the three surfaces |2](b)]. J^{z) for fiuorotalc is nearly identical to that for pyrophyl- 
lite, and that for talc features an additional local minimum around z = 5 A corresponding 
to water molecules above the binding site. However, the J^{z) curves at saturation are qual- 
itatively different from those at low RH [see[2](c-d)] For all three clays, and especially so for 
talc, the depth of the minimum at saturation is smaller than that at low RH, suggesting a 
weakening of adhesive interactions at saturation. 
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FIG. 3: J^{z) for adsorbing a single water molecule on the talc surface, compared to that 
for a molecule in the dimer and a molecule at saturation. 

To explore the competition between adhesive and cohesive interactions in talc, in [3| we 
compare J^{z) for an individually adsorbed water, with that for water in a dimer, and that 
for water at saturation. As shown in [3| the J-{z) for the dimer displays two minima. The 
minimum corresponding to the molecule inside the cavity is shifted to slightly larger values 
compared to the minimum in the J-{z) for a single water. In addition, the depth of the 
minimum is smaller, and is comparable to that for a single water on the more hydrophobic 
fluorotalc surface [2](c) and [3]. In other words, the presence of the second water weakens the 
adhesive surface-water interactions, which have to compete with the cohesive interactions 
between the waters. As the dimer is less tightly bound to the surface than a single water, 
it is easier for the water to escape the cavity in the presence of a second molecule. The 
dimer is in fact more mobile on the talc surface than isolated water molecules (not shown), 
confirming that the interaction of the surface with the dimer is weaker than with individual 
molecules. Finally, at saturation, cohesive interactions prevail, and water no longer occupies 
the binding site cavity as evidenced by the lack of a minimum in J^{z) for SA < z < 5A. 

F. Modified Talc Surfaces 

While the H-bonding between binding sites on the talc surface and water leads to an 
interesting transition from hydrophilic at low RH to hydrophobic at high RH, the binding 
sites interact weakly with water in fiuorotalc and pyrophyllite, which display hydrophobic 
behavior for all RH. To investigate the eflFect of the binding strength on the hydrophobicity of 
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FIG. 4: (a) J^{z) for a single water on various talc surfaces modified to span a range of 
A /iads- values, (b) The corresponding J^{z) curves at saturation, (c) The relative stability 
of water in the binding site compared to that in bulk, — A/Xgite, and the barrier to escape 
the binding site, A/Xbarrier, as a function of the binding strength, A/Xads- The dashed 
vertical line corresponds to /igat, the chemical potential at saturation. 



the surface, following Giovambattista et a/.^, we construct a series of modified talc surfaces. 
The only force field parameters that are changed are the charges on the oxygen (from 
qo = —0.95 to qo — 5q) and the hydrogen (from = 0.425 to + ^o) of the hydroxyl 
group. We study modified talc surfaces for 5q ranging from -0.425 which corresponds to a 
non-polar binding site similar to that in fiuorotalc, to +0.6 which corresponds to an ion-pair. 
5g = is the talc surface, by definition. 

In[4](a), we show J^{z) for an isolated water molecule on the modified talc surfaces. As 
the polarity of the -OH bond is increased, the magnitude of A/Xads also increases, providing 
us with surfaces that display a wide range of binding strengths. J^{z) at saturation, shown 
in[4](b) for these surfaces is particularly interesting. For weakly adhesive surfaces (—0.425 < 
5q < 0.1), there is only one stable basin at z 6.5 A, corresponding to molecules outside the 
binding site cavity. For stronger adhesion (larger a second basin develops at z 3.5 A 
and is separated from the first basin by a barrier. 

[4](c) shows the depth of this minimum relative to bulk, A/Xgite, as a function of A/Xads- As 
the surface becomes more adhesive, more waters occupy the binding site and the depth of 
this minimum increases. When adhesive interactions are large enough to overcome cohesive 
interactions, i.e., when — A/iads becomes larger than the chemical potential at saturation, 
— /Xgat (for 6q ^ 0.4), every binding site is occupied by a water molecule, resulting in a 
plateau in A/Xgite- 
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FIG. 5: (a) Schematic showing the surface coverage, B, over a wide range of relative 
humidities (RH = P/Psat ^ exp[/3(/x — /Xgat)]) and adhesive interaction strengths (A/Xads)- 
(b) Effect of A/Xads on surface hydrophobicity quantified by cos 9. The dashed vertical line 
corresponds to //gat- Snapshots indicating typical configurations of water molecules (red 
and white) on modified talc surfaces (blue) are also shown. As the adhesive interactions 
(A/Xads) overcome the cohesive interactions (/i), there is a transition from a dry surface 
[snapshots (i) and (iii)] to one covered with a monolayer of water [snapshots (ii) and (iv)]. 



However, the height A/Xbarrier of the barrier to escape the cavity, also shown in [4](c), 
continues to increase approximately linearly with the binding strength. Thus, for surfaces 
with strong binding, A/Xbarrier is large, and the exchange of molecules between the cavities 
and the liquid is expected to be very slow, with possible consequences on the extent of 
stick/slip at such surfaces in the presence of a hydrodynamic fiow. 



G. Tuning cohesion/adhesion via RH/A/Xads 

Collectively our results paint a comprehensive picture of how the experimentally mea- 
surable quantities, the surface coverage B, and the contact angle 0, respond to changes in 
relative humidity (or water chemical potential) , and on the strength of the adhesive surface- 
water interactions. The surface coverage B, is defined as the fraction of binding sites occupied 
by water molecules, and its dependence on RH and A/Xads is shown schematically in[5](a). 
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At low RH (= P/ Pqqx)) the competition between the adhesive interactions and the entropy 
of being in the vapor determines the surface coverage, B. At very low RH , there are no 
interactions between adsorbed waters and B can be approximated as : 

B ^ P/P'' = 0.1(P/Psat)e"^^^^^^"^•^ (2) 

where the second part of the equation is obtained by substituting for P* using Equation (1), 
and using appropriate values of the constants that depend on the surface geometry, amax 
and 5, and those that depend on thermodynamic conditions, T and Pgat- 

For surfaces with small adhesive interactions, i.e., — A/Xads < 5 kcal/mol (or — /5A/Xads < 
8.3), the coverage remains small (B < 0.1) even at saturation [snapshot (i) in [5]. Thus, 
no appreciable interactions between waters are expected over the entire range of RH- values. 
Both pyrophyllite and fluorotalc fall in this regime. 

Since B increases exponentially with /3A/Xads5 for values of — A/iads > 5 kcal/mol, there 
can be substantial coverage even at modest RH [snapshot (ii) in [5]. Equation (2) is then 
valid only for small RH- values for which the predicted B- values are small. Talc lies in this 
regime. 

For larger RH values, there are appreciable interactions between the waters, and it is the 
competition between adhesive and cohesive interactions that determines surface properties. 
For surfaces such as talc, for which — A/Xads < — Msat, cohesion prevails at saturation, and 
the adsorbed waters bead up into a droplet, while the rest of the binding sites on the surface 
are devoid of waters [snapshot (iii) in [5]. Thus, the interesting crossover from hydrophobic 
to hydrophilic behavior in talc is a result of its adhesive interactions being strong enough 
to overcome vapor phase entropy at low RH, but not strong enough to overcome cohesive 
interactions at saturation. In this regime, with increasing polarity of the binding site, the 
surface gradually shifts from hydrophobic to hydrophilic, and cos^ increases approximately 
linearly as shown in [5](b) . 

Finally, for surfaces with even larger values of — A/iads that are greater than — /Xgat, adhe- 
sion dominates. . Surprisingly, water does not fully wet the surface at saturation. Instead, 
all binding sites are occupied by water molecules and only this first layer of water wets the 
surface. This water is strongly bound to the surface and the microstructure of the surface 
dictates the relative positions of the waters. In the present case, the arrangement of waters 
on the surface is not commensurate with the hydrogen bonding network of water, so that 
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water beads up on the monolayer [snapshot (iv) in [5) . For the modified talc surfaces with 
— A/Xads > — Msat, the surface has a strongly adsorbed water monolayer with a droplet on it 
that makes a contact angle of about 50°. 

Similar behavior was reported by Ohler et al for titanium dioxide surfaces, with droplet 
contact angles of 32 — 34° on top of roughly two monolayers of watei^^. However, other 
simulation studies investigating the effects of surface polarity on hydrophobicity'^^'^, do not 
observe a plateau with non-zero contact angle at large polarities, seen in our results (5](b)]. 
Our modified talc surfaces are different from these previous studies in that the variation in 
polarity was achieved by changing the charges on atoms in recessed binding sites, while the 
remaining surface atoms remained the same. In contrast, in ref.'^, the surface was modified 
by changing dipoles that protrude from the surface, while leaving the remaining surface 
atoms unchanged; whereas in ref.'^, the charges on all atoms in the top two layers of an 
FCC crystal (111 facet) were changed to tune the polarity. Thus, our results indicate that 
the microstructure of the surface is important in determining the effect of polarity on its 
wetting properties. 

In contrast to the wetting properties of the model FCC surfaces used in ref.'^, experi- 
mental measurements indicate that the FCC crystals of platinum (Pt), palladium (Pd), and 
gold (Au) are hydrophobic. Kimmel et al observed a hydrophobic water monolayer on 
both Pt(lll) and Pd(lll) crystal^^ZlSI Similarly, water has been shown to bead up on Au 
surfaces^^ with a contact angle of 100° and Au surfaces have also been shown to adsorb, and 
facilitate the unfolding of proteins^^; behavior that is typically associated with hydrophobic 
surfaces^^. We speculate that the hydrophobicity of these metal surfaces arises from the 
presence of a monolayer of water, which binds strongly to the surface in a geometry that 
inhibits hydrogen bonding to the subsequent liquid water molecules. 

Our results also have implications on the wetting properties of nanoporous silicates such 
as hydrophobic zeolite iJ^^^ti^ and metal-organic frameworks^^. These hydrophobic pores are 
thought to be devoid of water at ambient conditions, with water intrusion into the pores 
occurring only at sufficiently high water pressures. Our results suggest that in the presence 
of strong binding sites, these nanoporous materials may contain strongly adsorbed water 
molecules, even at lower pressures. If the resulting water-covered surface is hydrophobic, no 
further filling of the pores (analogous to wetting for planar surfaces) would be observed at 
ambient pressures, and intrusion would occur only at higher pressures. 
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Appendix A: Unit cells 

The unit cells used for the simulation of talc and pyrophyllite are reported in |6] and [7| 
For ffuorotalc, the oxygen of the hydroxyl group is replaced by a ffuorine atom and the 
hydrogen is removed. The unit cell of pyrophyllite, a dioctahedral smectite, has dimensions 
along the surface of 5.18 x 8.97 A^, as known from X-ray diffractioiF^. The unit cell of 
ffuorotalc is not known exactly; we used the one determined by X-ray diffraction on synthetic 
ffuorohectorite^^, which differs from ffuorotalc only by substitution of some magnesium by 
lithium in the octahedral layer, resulting in a permanent negative charge compensated by 
sodium counterions. The unit cell has dimensions 5.24 x 9.09 along the surface. For talc 
we used the same structure, replacing each ffuorine by a hydroxyl group with a bond length 
of 1 A, oriented perpendicular to the surface. 

Appendix B: Comparison of force fields 

In the present work, we used the CLAYFF force field to describe the clay surfaces and their 
interactions with water molecules. To justify this choice, here we compare the predictions of 
another commonly used force field, and those of CLAYFF, with experimental results. This 
force field was originally developed by Skipper et al^^ and adapted by Smith et al^^ for its 
use in conjunction with the SPC/E water model. 

To investigate the talc surface at low RH, in[8](a), we show the J^{z) obtained using the 
Skipper/Smith force field and compare it with that obtained using the CLAYFF force field. 
Also shown is the experimental estimate discussed in the main text, indicating that the 
Skipper/Smith force field overestimates the binding or adsorption free energy. 
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To investigate the hydrophobicity of talc surfaces at saturation, obtained using the two 



Skipper/Smith talc surface indicates that it is harder to empty the observation volume 
close to the surface than in bulk water. This is also consistent with the observed complete 
wetting of the talc surface by a droplet, indicating a contact angle of ^ = 0°. Such a complete 
wetting is however in contradiction with the experimental contact angle of 80 — 85°. We thus 
conclude that the Skipper/Smith force field significantly overestimates talc- water adhesive 
interactions, both at low RH and at saturation. 

Another force field used to model dioctahedral clays and their interaction with organic 
cations was proposed by Heinz et al.^^. This model was not extended to triocahedral clays 
such as talc, and the behavior of water at clay surfaces modeled with this force field has 
not been reported. We nevertheless simulated water droplets on the surface of pyrophyllite 
using this force field. The resulting contact angle (125°) was larger than that measured 
experimentally (80 — 85°), suggesting that this force fields results in surfaces that are too 
hydrophobic. 

Finally, while we find that CLAYFF is the best available force field to date, to simulate 
water at the surface of uncharged clay minerals, the present work suggests that it is too 
hydrophobic. Thus we find that there is room for improvement to describe the clay- water 
interaction, in agreement with the findings of a recent study comparing molecular simula- 
tions with X-ray and neutron diffraction experiments on a charged smectite^^. The insights 
gained during the present study of neutral clays, which are more sensitive to the clay- water 
interactions, could also be helpful in the design of an improved force field. Such design 
requires a subtle balance between diflFerent interactions which is generally not achieved by 
tuning only one parameter. With this caveat in mind, we note that a slightly more polar 
hydroxyl group might be relevant, as the modified talc surface with 6q = 0.1 seems to agree 
quite well with experimentally measured A/Xads and cos 9 values for talc. 
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FIG. 6: Atomic coordinates in the talc unit cell. Subscripts for oxygen differentiate 
tetrahedral (Td), bridging (B) and octahedral (Oh) atoms. 
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FIG. 7: Atomic coordinates in the pyrophyllite unit cell. Subscripts for oxygen 
differentiate bridging (B), tetrahedral (Td) and octahedral (Oh) atoms. 
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FIG. 8: (a) J^{z) for the adsorption of an isolated water molecule on talc simulated using 
the CLAYFF and Skipper/ Smith force fields. The arrow indicates the experimental value 
of the minimum, estimated by fitting the adsorption isotherm of Michot et al^ to a 
Langmuir model in the very low RH regime (see text), (b) Pv{N) for the talc surface, 
using the CLAYFF and Skipper/Smith force fields. 
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